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Abstract 

Full manifold of the complex Bloch-Floquet eigenfunctions is inves- 
tigated for the ground level of the purely magnetic 2D Pauli operators 
(equal to zero because of supersymmetry). Deep connection of it with 
the 2D analog of the " Burgers Nonlinear Hierarchy" plays fundamen- 
tal role here. Everything is completely calculated for the broad class of 
Algebro-Geometric operators found in this work for this case. For the 
case of nonzero flux the ground states were found by Aharonov-Casher 

(1979) for the rapidly decreasing fields, and by Dubrovin-Novikov 

(1980) for the periodic fields. No Algebro-Geometric operators where 
known in the case of nonzero flux. For genus g = 1 we found periodic 
operators with zero flux, singular magnetic fields and Bohm-Aharonov 
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phenomenon. Our arguments imply that the delta-term really 
does not affect seriously the spectrum nearby of the ground 
state. For g > 1 our theory requires to use only algebraic curves with 
selected point leading to the solutions elliptic in the variable x for KdV 
and KP in order to get periodic magnetic fields. The algebro-geometric 
case of genus zero leads, in particular, to the slowly decreasing lump- 
like magnetic fields with especially interesting variety of ground states 
in the Hilbert Space £ 2 (M 2 )- 

Introduction. Magnetic Pauli Operator and 
factorizable Schrodinger Operators 

A nonrelativistic 2D Pauli operator for the charged particles with spin=l/2 
moving in electric and magnetic fields E a = d a U, A a (under the Lorentz 
gauge condition) has a form (see [T], let e = 1, m — 1/2, we neglect the 
universal constants c, H whose values are inessential here ) 

2 

a=l,2 a=l 

a a are the Pauli matrices 

ai = ( i o ) ' a2 = ( < o* ) ' a3 = ( -1 ) ' a ° = ( 1 ) • 

Let U = (i.e. electric field is equal to zero). The operator L p is reduced 
to the direct sum of two scalar Schrodinger operators. They are written in 
the "factorized" form (see S] ) : 

L p = QQ+ © Q+Q = L + ®L_ (2) 

(we neglect all constants unimportant for our goals). Here Q = d + A, 
Q+ = -(B- A) ± -A = A&, dA + OA = 2B is magnetic field, d = d x - id y , 

8 = d x + id y , dd = a = d 2 x + d 2 y . 

Magnetic field B here is perpendicular to the oriented plane (x, y). There- 
fore, the magnetic field has a sign. For the rapidly decreasing field B at 
x 2 + y 2 — > oo, the magnetic flux {B} is finite by definition. The Operator L p 
is nonnegative. Therefore, the ground state energy is equal to zero Eq = 



2 



or positive Eq > 0. Let \{B}\ > 1 (in natural quantum units). Then the 
ground state subspace in the Hilbert Space is a linear space of the dimension 
[{£>}] = m (see [2]). For the periodic case and integer flux {B} G Z, the lin- 
ear subspace of ground states is infinite-dimensional. It is isomorphic to the 
Landau level in a homogeneous magnetic field (see [31 E]). The higher levels 
are separated from the ground state by the nonzero gap A# > 0. According 
to the literature of the late 80-s, this operator admits a "supersymmetry" 
transformation P : L + ->■ L_ 0, P 2 = 0, P : ^ -> Q + \fJ for {B} > 0, 
iff G L + ; All positive energy levels e > are degenerate since (\I/,P\I/) both 
are eigenfunctions. For the zero energy we have P^f = 0, Eq = 0, if function 
iff belongs to the Hilbert space £ 2 (M 2 )- Using the ancient language of XIX 
(or even XVIII) century, there exists a "Laplace transformation" of the 
2D second order scalar operators (see [TT] ) 

L = (d x + A){d y + D) + U, U = e f 

L^L = e f (d y + D)e~ s (d x + A) + U (3) 

The equality L^f = implies = 0. We have here L = QiQ 2 + U. 

In the selfadjoint elliptic case, studied in detail in [TT] from the point of 
view of spectral properties of operators, we have 

L = QQ + + U, U = e f 

Q = (d + A), Q + = -(d-A). (4) 

For the "purely factorizable" operators U = const, we have L = Q + Q + U 
(let U = 0). The Laplace Transformation 4? = coincides with the 
"Supersymmetry" with P = Q + in the sector L + and P = in the second 
sector L~ . It acts on the whole spectrum = e^, = e^>. The Ground 
States of L p are all "instantons" , i.e. they satisfy to the equation Q + ty = 
if G £2(JBL 2 ). Let zero mode L^f! = does not belong to the Hilbert space 
but "belongs to the spectrum". It simply means that its growth rate is 
slower than some polynomial for x 2 + y 2 — > oo,. The instanton argument 
disappears. In the last case the point e = is the bottom of continuous 
spectrum. Even if no "instanton" type solutions ^ of that kind exist, it 
is impossible to conclude immediately that the true ground state for the 
operator L is positive Eq > (but it is highly probable). 
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In the case of nonzero flux {B} ^ the ground energy of the Pauli 
operator is equal to zero e = 0. For L p = L + © L_ it is realized inside of 
the sector L = L + (if {B} > 0) or L = L_ (if {B} < 0) (see S S]), in 
the rapidly decreasing and periodic case (see also [5] , where other functional 
classes of magnetic fields were considered). The rest of the spectrum is 
twofold and separated from zero by a positive gap A B . Interesting classes of 
the "factorizable" operators L, having one more infinitely degenerate level 
except of the ground one, were found in [II] . These works have a "soliton" 
origin. Let us point out that the connection between Laplace transformations 
and 2D Toda chain found in the soliton theory, was in fact known in the XIX 
century to Darboux and his school. However, all calculations at the end of 
XIX — beginning of XX centuries were purely formal, and the elliptic case 
was completely missing. 

In the present paper we investigate the algebrogeometric case. For the 
smooth periodic operators (i.e vector-potentials are periodic) we have mag- 
netic flux equal to zero, but for the degenerate soliton-type case it might 
be not so. In our case the whole complex variety T of the Bloch-Floquet zero 
level eigenfunctions L\l/ = appears ("The Complex Fermi Curve"). This 
operator is called "algebro-geometric" if genus is finite. 

A Purely Factorizable Reduction of the self-adjoint Schrodinger 
operator L = — (d + A)(d — A) is studied here from the point of 
view of algebro-geometric operators. It was recently found by the 
present authors [T2] using the 2D Soliton- Type Completely Integrable Sys- 
tem. The operator L = d x d y + Gd y + 5 is by definition hyperbolic in the 
work [12]. Its reduction 5 = leads to very interesting "2D Burgers Hier- 
archy" which is linearizable similar to the classical Burgers Equation. The 
spectral meaning of this linearization is revealed in [T2]. For the nonreduced 
system L t = (LH — HL) + fL the second operator H = A + Fd y + A 
may have an interesting spectral theory in the stationary "finite-gap" case 
LH — HL = —f(x, y)L. The operators L and H form "The Algebrogeometric 
Pair of PDE's Commmuting Relative to the Level of Operator L = 0" accord- 
ing to the terminology used by Krichever and Novikov in the late 1970s-early 
1980s. The operator H is elliptic here. Its study is the second main goal of 
[12]. It is easy to make H smooth, periodic and real. However, we failed to 
find nontrivial self-adjoint operators H within this approach. So the Con- 
jecture is formulated in [12] : For the smooth periodic self-adjoint 2D 
Schrodinger operator in R 2 the Full Complex Manifold of Bloch- 
Floquet Eigenfunctions might contain Complex Algebraic Subman- 
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ifolds only belonging to one energy level (except some trivial cases 
which essentially can be reduced to one-dimensional equation). 

1 Algebro- Geometric self- adjoint factorizable 
operators. The inverse spectral data 

As it was demonstrated above, the purely magnetic 2D Pauli Operator L p = 
L + ©L_ is a direct sum of two "Strongly Factorizable" Schrodinger operators 

L + = QQ + , L_ = Q + Q 

Q = {d + A), Q + = -{d-A). 

Following [IB] , let us recall what is the algebro-geometric operator with pe- 
riodic coefficients A, U 

L = -(d + A)(5 + D) + U, 

where A, D, U are periodic in x, y. 

Let us describe "The Inverse Spectral Data" for the Operator L: we take 
nonsingular Riemann surface V of genus g > 0, two marked "infinite" points 
ooi, oo 2 with local parameters fc /-1 (ooi) = 0, fc //-1 (oo 2 ) = 0, and set of g 
points ( a "divisor" of degree g) V = (Pi, . . . ,V g ). We write it as a formal 
sum T> = V\ + . . . + V g . In the work [13] the "Two-point Baker- Akhiezer 
function" was introduced, *$?(V, x, y), V G T, with the following properties: 

a) It is meromorphic in the variable V outside of infinities; 

b) It has following asymptotic near infinities oo 1; 002 : 

oo, : * = c(x, y y (l + O (±)) , i ( oo,) = 

°* = *=«■"• + ° Of)) ■ fR^ 

c) It has poles of the first order in the points of divisor T> which are inde- 
pendent on the space variables x, y. 

d) * = 1 at x = 0, y = 0. 
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Such function satisfies to the equation 

L^ = 0, L = dd -2(\nc) z B+U(x,y). 

Our function \1/ is a Bloch-Floquet eigenfunction if c(x,y), U(x,y) are pe- 
riodic. In general, they are quasiperiodic but for a dense set of Riemann 
surfaces T these functions are periodic. It is possible to write down c, 
U through the 0-functions in a standard way of the periodic soliton theory 
(see [13] and surveys [IB], [T7]). A Purely Potential Reduction c = 1 was 
found for this data in [13], [IS]- The Self-adjoint reduction in the presence 
of magnetic field B = — A(lnc)/2 = — <9<9(lnc)/2 ^ was found in [18J: An 
Antiholomorphic Involution a : Y — > T, a 2 = 1, cr(oo!) = oo 2 must exist such 
that 

a{k') = -F, a{V) + V ~ (K) + oo! + oo 2 , (5) 

where (K) is divisor of zeros and poles of holomorphic 1-forms, and the 
symbol ~ means "linear equivalence " of divisors, i.e. every divisor of zeros 
and poles of meromorphic function is equal to zero. These conditions are 
sufficient and (probably) necessary, but no rigorous proof of necessity was 
obtained in the literature yet. 

Let us describe "The inverse spectral data" for our " Factorizable 
Self-adjoint Reduction" . It is result of the present work. 

Riemann surface T is degenerate 

r = r' u r", 001 e V, oo 2 e r", 

and the intersection V R T" is a set of k + 1 points Q' Q ,. . . ,Q' k . 




vanishing cycles vanishing cycles 

Fig. l.£ = Fig. 2. g = 1 
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An antiholomorphic involution should exist a : V — > T", T" — > V, which 
permutes points ooi ->■ oo 2 , oo 2 ->■ ooi, (O/ , • • • , Q' k ) ->■ (<5L • • • , <5',)- Let 

g = genus T' = genus r". Let us specify g + /c points X>' = (V[, . . . , V' g+k ) on 
T' and g points V" = (V", . . . , Vg) on r" satisfying to the linear equivalence: 

a(V + T>") + V + D" ~ (if) + ooi + oo 2 on T, (6) 

where (K) = (K') + (if") is a divisor of 1-form 00 with conditions on residues: 

w = J ( r') 
w = u/' (r") 

-Res w' = Res w", crfi?) = Q. 

Qj Rj=a(Qj) 

Let us reformulate it in terms of one curve V with local parameter 1/k', 
set of points Q ,. . . ,Q n and a(Qj) = Q a (j)- 
Find \l/-function on V such that: 

1) It has poles of first order in points V, (g + k points). 

2) It has asymptotic near 00 1: 

* = c(x,y)e k ' s (l + O 

3) \P| = (pj(z) holomorphic in z. 

Q'j 

(k + 1) points Qq = Q' io ,. . . ,Q" k = Q\ are fixed on the surface T'. Let us 
construct an antiholomorphic (one-point) Baker-Akhiezer function tp with 
properties: <f(V, z) is antiholomorphic in V G V; it has a fixed divisor of 
first order poles cr(V") and asymptotic 

<p~e-* z (l + o(jj;yj, <p=l x = 0,y = 0. 

It is necessary to satisfy the condition 

a{V") +V = (Q) + (K')t> + 001 

with restriction on residues (above) in the points Qj ~ Qa(j) f° r the form 
defining (K 1 ). 
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Let 

<p j (z)=<p(Cy j ,z), Q" 3 =Q' a{]) (7) 

We described everything in terms of one nonsingular surface V because Y" = 
nil"). 

Under these conditions we prove that our function ^f(V, z, z) on the sur- 
face V 3 V satisfies to the equation below in the space variables z, z: 

L + ^ = 0, L + = (d + A)d, 

where A = —d\nc. There is a constant q^O such that the function (etc) is 
real. It generates symmetric operator by the following rule: let 

c = e 2 *, L + = -e"*(0 + A)de* = 

= -(d + A)(d - A) = QQ + , A = A/2 (8) 
* ^ e -*^ = J-^, i ^ 

The operator L + is nonsingular, if there exists a constant a such that ac(x, y) > 
0. Let 5o — be a point of the spectrum for L + . A natural candidate for the 
ground state is a function \j \fc satisfying to the equation 

Q + (1/Vh) =(-B + A)e-* = 0. (9) 

As (d + A)e q> = 0, the function e* is also a candidate for the ground state for 
the operator L_ = Q + Q. However it is true only if these functions belong to 
the spectrum in C 2 (M. 2 ). Let us formulate sufficient conditions for that: 

a) The coefficients are periodic (quasiperiodic) and nonsingular ac > 0. Both 
functions \/ac, \J\jac are positive and belong to the spectrum. We know 
that A = — (<9mc)/2, B = —(A /nc)/2, and the magnetic flux is equal to 
zero: 

= / / Bdz A dz = 0, where □ is an elementary cell. 



It is true because B = —(dd\nc)/2, Bdz Adz = —& Z zdzAdz = —d(Q z dz) 
1/2 • d(Adz), so our flux is an integral of the exact form. 
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b) Let us consider an exponential case now. 

Ground States of the Nonrelativistic Pauli operator: 

Let L = L + © L_, L + = QQ+, L_ = Q+Q, where Q = (d + A), Q+ = 
(-d + A), and A = -2$ 2 , A = -2$ 2 , where $ is real. 

In the whole class c — > e w c = d, where W = ax+(3y for real a, (3, we have 
the same magnetic field B = -A$ = -(Alnc)/2 = -(A In (ce w ))/2. 
The operators L + and L' + are unitary equivalent. Indeed, we have: 

L ± = -(d x -t%) 2 -(d y + i$ x ) 2 ±A$ (10) 
L' ± = -{d x - t% - t/3/2) 2 - (d y + i$ x + ia/2) 2 ± A$ (11) 

The Unitary Gauge Transformation 

l]/ ^ r gi(-/3a;+a J /)/2 _ ^/ ^ 2 ) 

L+ ->■ L' + , L_ ->■ L'_ (13) 

realizes this equivalence. If any one of the functions \fd or 1/ \fd is 
bounded, then e = 0, and this function serves the spectrum in £ 2 (M 2 )- 
We have = e'^-^y/d, for L+, ^_ = e'^-^l/Vd for L_. 
Thus we constructed as many different ground state vectors Eq = for L + 
or L_ as there are bounded functions {e^c} or {e~ w c~ 1 } in this class. 

Example 1. Consider c = 1 + e y . It does not depend on x. Here c _1 
is bounded, and c is unbounded. Let d a g = e ax+l3y c. To have bounded 
function 1/c' in the class d = e w c,W = ax + j3y, we need to satisfy 
the conditions: a = 0, — 1 < f3 < 0. So we obtain continuum of 
eigenfunctions for the ground energy level e = parametrized by 
the index f3: 

V^e^-l/Vd, (5e [-1,0] 

Here Co t o — c, — 1/y/c. The unbounded functions (\/c, 1/y/c) 
satisfy to the equation L + (y/c) = 0, L_(l/y/c) = but do not 
belong to the spectrum. 

So we are coming to the following conclusion: 

For the purely exponential generating function c > 0, the level Eq = 
is the lowest point of the spectrum if we can find in the class {e ax+ ^ y c} 
for real a, f3 a bounded function d = e ax+(3y c or 1/c' = e~ ax ~P y /c, L p = 
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L + © L_ . But it is always true for all k > where k + 1 is the number of 
intersection points if all coefficients positive. 

In the smooth periodic case we know that the smooth periodic functions 
y/c and l/y/c both are the ground states if c is positive. They are peri- 
odic, and the zero energy level is always a twice degenerate point of the 
spectrum (not like in the case on nonzero flux). 

In the next paragraphs we present calculations for the genuses g = 0,1. 

2 The Algebogeometric self-adjoint factor- 
izable operators 

2.1 Solutions of genus g=0. 

As we can see below, all algebrogeometric purely magnetic Pauli operators 
with Complex Fermi Surface of genus zero correspond to the functions c 
of the form 



Here constants, and Wj = ajz + bjZ are the linear forms with 

constant coefficients. In general, all coefficients here are complex. How- 
ever, for the selection of physically meaningful self-adjoint operators and 
real magnetic fields B = — (Alnc)/2 we are going to formulate proper 
restrictions below. These "Degenerate Algebrogeometric Purely 
Magnetic Pauli Operators" are the n Natural Analogs of the 
"Multisoliton Potentials" for the ID Schrodinger Operators in 
the case of KdV. It deserves to point out that our generating functions 
c are linear combinations of the elementary exponents with constant co- 
efficients. Indeed, the magnetic field B is equal to — l/2(Alnc), i.e. it 
is a strongly nonlinear object. Such linear behavior of the quan- 
tity c reflects the main property of the "2D Burgers Hierarhy" 
discovered in the work |12j . For KdV the Multisoliton functions are 
also constructed as the second logarithmic derivatives of something which 
is indeed a nonlinear expression like some determinant made out of the 
one-soliton functions. In spite of linearity of c, the Spectral Theory 
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is quite nontrivial for these 2D Purely Magnetic Pauli Analogs 
of the Multisoliton Operators. 

In the case of genus zero our \I/-function is written in the form (k — k' e 

n 

c = u ,* = e'« , U ° kn + '" + Un r , V = (a 1 ,...,a n ) (11) 
(k - ai) . . . [k - a n ) 

(Here n + 1 is the number of intersection points) 

r = r u r", r' = r" = s 2 

with local parameters k = k'(T') p = k"(T"); points 001,002 have the 
form k = oo(ooi), p = 00(002). We have <p(z,p) = e pz . The intersection 
points are fc , . . . , k n for V and po, . . . ,p n for T". They lead to the set of 
equations for ^ in these points: 



V\k=k s = e p > z ,j = 0,...,n. 
So our solution has a form: 

A( n_1 ) 



3=0 



(12) 



(13) 



as it follows from the system of linear equations 

uokS + --- + u n = (ko-a 1 )...(k - a n )e PoZ - k ° z 
u k% H h u n = (k n - ai) . . . (k n - a n )e p " z ~ knZ . 



(14) 



Here 



A (n) 



kg ... 1 



k r : 



^i<j{ki kj) 



and A^ n ^ are similar Vandermonde determinants with the set of generat- 
ing numbers (ko, . . . , kj, . . . , k n ), where kj is erased, 9j = (kj — ai) . . . (kj — 
On)] u = c. 

The quantity c(x, y) = u Q is determined by the field B up to the trans- 
formation c — » ae w c = d , where a = const, W = jz + 5z because 
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B = -(Alnc)/2 and -(Alnc')/2 = -(Alnc)/2, A = dd. Therefore we 
have exactly n unknown coefficients in the formula (15): 

n 

j = 0,...,n, c = Y,^ Wq{z '"\ (15) 

q=0 

For the given kj,pj all coefficients determined up to the common 

multiplier. 

For the differentials below the conditions on residues should also be sat- 
isfied: 

^ _ (k - oi) ... (k - a n )dk _ s(p + ai) . . . (p + a n )dp 
1 (k — k )...(k — k n ) 2 (p-po)...(p- Pn) 

where s is a constant. 

ReSfc.fii + Res Pj f2 2 = 0. 

Choosing appropriate divisors V = ( ), we obtain all complex 

coefficients Kj G C. 

We need to classify such divisors V and linear forms Wj = pjz — kjz that 
c is real and positive in the equivalence class c — > ac = d, a = const. 

For the reality of c, (x,y) & R, z = x + iy,z = x — iy, our expression must 
consists of the following terms: 

1. The "exponential type" term for some index j: 

pj = -kj, Wj = pjZ + pjz, 

Kj is real and Kje Wj is also real (a purely real exponent) 

2. The "mixed type" term for the pair of indices (j, I): 

Pi = -kj,k t = -pj, 

R,j = Ki, Kje Wj + Kie Wl = Kje Wj + Rie Wl is real. We assume that all points 
k q , j 7^ q and pj ^ p q , j ^ q are distinct. 

For the case 1: We obtain terms like real exponent Kje^ ajX+ ^ jy \ where 
Pj = otj + i(3j, kj = —pj, Kj — Kj = 0. 
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For the case 2: We obtain terms of the form 

e w ™ (k'j cos W Id - k] sin W Id ), Kj = k'j + 

Wj = Wrj + iWjj = [(a, - -Yj)x - % - 5 3 )y] + - Sj)x + (a, + ^)y], 
where pj = —k q = aij + ij3j, kj = —p q = jj + iSj. 
For Wjj = 0, we have the case 1: pj = —k. 

3. The "purely trigonometric type" appears as another special subcase 
of the case 2 if W R j = or atj = Tj,/^ = —5j, i.e. kj = pj, Wjj = 

(/3jX + ajy),e w i = 

In all these cases c is real. The mixed case 2 leads to the zeros of c and 
singularities of magnetic field if they are not " blocked" by other stronger 
terms. 

c = J2 Kje Wj{x ' y \ c -)• ne w c = c', 
j 

where all Kj are real. Let c = c + + c~. Here Kj > for j e I, Kj < 

jell. 

Consider at first the case c = c + , i.e. Kj > for all j. So we have c > 0. 
As one can see (below), the magnetic field B = — (Alnc)/2 is bounded 
in M 2 . In the class {ne w c} both \fd ', ^= never can be bounded. Either 
they both exponentially increase along some directions at x 2 + — >■ oo 
or one of them (i.e. (l/y/d)) is bounded. In the last case the pair {c, W} 
or simply a function e w c = d defines the ground state. As we can see 
below, such functions d form a domain inside of the convex polygon T in 
R 2 which is always nonempty. This domain is completely determined by 
the set of linear forms Wj in the class {VFj}: it is a convex hull of the set 
of points (ccj, Pj) G R 2 , Wj = atjX + fay (see below). 

2. Let us consider the purely trigonometric case. Here we have the 
cases of odd and even numbers of intersection points n + 1. They are 
drastically different. 

a) The number of intersection points n + 1 is even. 

b) The number of intersection points n + 1 is odd. 
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a) c= /tj cos Wij + k" sin Wij, 

3=0 

where Kj = k'j + iK",Wj = iWij = —kjz — kjz,kj = otj + i(3j. Here 
function c always has zeros. 



b) c = 1 + cos Wjj + /c" sin PF 7 



For the appropriate constants k', k" G R we have c > and magnetic 
field £> = — (Amc)/2 is smooth, periodic and has zero flux through the 
elementary cell of periodic lattice (or the quasiperiodic mean value if c is 
quasiperiodic). It would be interesting to describe corresponding domains 
in the space of constants. The set of lines Wjj = ol 3 x + fyy, should pass 
through the integer vectors of the lattice in IR 2 . Otherwise, the fields are 
quasiperiodic. Our conclusion is that in the regular trigonometric case 
both functions y/c, \j\fc are periodic and positive; they are the 
ground states in both sectors of the Operator L p . In the 

general quasiperiodic case the situation is the same. 

A Curious Remark. There are "critical" values of constants k'- k"- such 
that c has isolated zeroes c = (repeated periodically). It is possible to 
choose parameters such that we have in this critical point an isotropic 
hessian d 2 c = ±a 2 (dx 2 + dy 2 ). Then the magnetic field has a 5-shaped 
singularity B = — (Alnc)/2 ~ 5(x — x ,y — y ). 

Example 2. Let n = 4. We demonstrate here a simplest nonsingular 
purely trigonometric (i.e degenerate algebrogeometric) operator, essen- 
tially dependent on both variables x, y: we write \l/-function in the form 

* = e {k 2 - a\){k 2 - a 2 ) ^=K«2>-«i>-° 2 )> 

and (f = e pz . 

Take the intersection points of V and T" in the form 0, ki, k 2 , —k±, —k 2 
for T' and 0,pi,p 2 , — Pi, —?2 for T". Let 

Pi = fci G R, p 2 = -k 2 = %K G iR. (16) 
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In this case the antiinvolution a : k — > —p is correctly defined on T. The 
differentials look like 

_ (g - aj)(k 2 - a\)dk 

1 (k 2 - kj)(k 2 - k 2 2 )k ' 

= s(p 2 -af)(p 2 -g|)rfp 

2 (p 2 - pI)(p 2 - pDp 

where s is some number. The condition on the residues 

Res ^i + Res ^2 = 0, Res±fc.fii + Res± Pj Q 2 = 0. 
must be valid. In the points of intersection we have 

(0) = 1, ^(h) = e piz , V(k 2 ) = e P2Z , ^(-h) = e- piz , V(-k 2 ) = e^'. 
So the equalities follow: 









2 2 

W4 = a 1 a 2 , 








u kf 


+ u\k\ + u 2 k\ 


+ u 3 h 


2 2 , / 7 2 

= — a 1 a 2 + (k x - 


-aim 


2 

- a 2 




uok 2 


+ u\k\ + u 2 k\ 


+ u 3 k 2 


2 2 , / ; 2 


-a\){k 2 2 


- «2 




u Q k\ 


- u\k\ + u 2 k\ 


- u 3 k x 


= —afa 2 , + (A; 2 — 


<*?)(*?- 


-a 2 ) 




u§k\ 


— u\k\ + u 2 k 2 


- u 3 k 2 


2 2 , / ; 2 

= -0^2 + {k 2 - 


a 2 ) (A; 2 - 


-a 2 ) 





Sum of the second equality with the fourth equality and of the third 
equality with the fifth one are written below: 

2u kf + 2u 2 k 2 = -2a\a\ + {k\ - a\){k\ - a 2 2 )2 cos(2hy), 

2u K 4 - 2u 2 K 2 = -2a\a\ + (K 2 + a\)(K 2 + a\)2 cos(2Kx). 

Take 

di G M, a 2 = ia G iR. 

Using (16) and taking s = —1, we can see that the conditions on residues 
of the differentials Qi and Q 2 are satisfied. We have 

2 2 

c = u = {I- A cos(2k iy ) - B cos(2Kx)) , 
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where 



(k 2 -a 2 )(k 2 + a 2 ) (K 2 + al)(K 2 -u- 



k 2 (k 2 + K 2 ) ' ~ K 2 (k 2 + K 2 ) 
For = 10, fci = 5, ai = 2, a = 1 we obtain 

546 „ 2574 „ „ 624 
5 = , A + 5 



3125' 3125' 625' 
So the magnetic field is smooth and periodic. 

Example 3. Consider now the Full Class of the Purely Exponen- 
tial Real Type functions C. We introduce below an important notion 
of "the Indicator of Growth" for the set {Wj} of all real exponents en- 
tering formula for c with positive coefficients. Suppose this function c 
and therefore yfc grows exponentially in all directions. This is a Stable 
Property. The function ^ has exponential decay in M 2 . Many other func- 
tions l/\fd in the same class d = e w c are such that l/y/d automatically 
have similar decay (for example, it is true for all "small" linear forms 
W = e(ax + by),e ->■ 0). 

Consider first the Unstable Case n — 1, i.e. with 2 intersec- 
tion points. Let us take absolutely typical example c = 1 + e y follow- 
ing notations in the Example 1 above. The polygon T numerating all 
bounded functions 1/d in the class d = e w c, coincides with a segment 
(3 G [— 1, 0], W — ax + (3y. It does not have inner points. We never have 
d in this class which has exponential growth in all directions. Magnetic 
field here depends on one variable. The same result is true for all 
cases with n + 1 = 2 where c = K\e Wl + n 2 e W2 . The Unstable cases 
for all n > 1 are given by the sets linear forms {Wj} such that 
all differences are proportional to one linear form with constant 
coefficients. Magnetic field here depends on one variable only. 
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R 2 = (a, 13) 
Wj = ctjX + fyy 





P 




-^=(0,0) 




a 




T 


— ^» 


-W 2 =(0-1)' 






Fig. 3 





Let us consider examples of the indicators of growth in this class. Choose 
W = e~ y l 2 . We have d = e y l 2 + c~ y / 2 . The indicator of growth has 
zeroes I w i(<p) = exactly in two points ip = 0, 7r. Put W = e^ x ~ y ^ 2 . We 

get c" = e^ x+y ^ /2 + c( x ~ y ^ 2 . For this case Iw"(f) = on the connected 
segment. However, the zero set of the indicator of growth is never empty 
for n — 1. It is not surprising because magnetic field always depends on 
one variable for n — 1. 

c=ey /2 + e-y /2 




Fig. 4 a. Fig. 4 b. 
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(the growth directions for Wj are shown by the rows; every ray from the 
zero point belongs to at least one sector of exponential growth provided 
by the 3 exponents entering c.) 

The Stable Cases start with n = 3 (like in Fig 4,5). Let, for example, 

c = e Wl + e W2 + e Wz = e x + e y + e~ x ~ v 



a) R 2 = (x, y) 



b) 
W 3 



l 2 = (a 7 (3) 




W x =x (o=l, p=0) 



|W 2 =y (o=0, p=l) 



"P -w 3 =(i,i) 



-\Hf=(-l,0) 



Fig. 5 a. 




In this example the function 4 = ~w- is bounded if and only if W E T 
where T is a triangle with vertices (—Wi, —W2, — W3) C R 2 with coordi- 
nates (a,/3). It belongs to the Hilbert Space (i.e. it is square integrable 
function on the x, y-plane R 2 ) if and only if d = e w c where W belongs 
to the interior of T. 

We describe below all nonsingular cases for the exponential gen- 
erating functions c: 

Definition. Let W = ax + f3y be a real linear form. We call the function 
on the circle Iw > 0, Iwiv) — max(acos(^ + (3 sirup, 0) "the Indicator 
of Growth" for the linear form W = ax + f3y. For the set of real 
linear forms {Wj} we call function I{ W .y((p) = max(7vi^) > the "the 

Indicator of Growth of this set". It is what some people call "a 
Tropical Sum". 
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I. Consider the case 

ci = e Wj Kj, Kj > or en = e Wj Kj, Kj < 0. 



Let us note that the indicators of growth are different for the 
different sets within the same class {Wj} and {Wj + W}. It is 
an invariant of the set of exponents entering the function c, not 
of magnetic field. It does not depend also on the coefficients Kj 
entering the function c. 

There are following possibilities: 

1. I{ Wj y(tp) > for all ip G S 1 (see Fig. 6 a). 

2. I{ W .y((p) = on the connected closed segment on the circle (p G S 1 
(see Fig. 6 b). 

3. I{Wj}{v) = in the isolated points ip = ipi, . . . ,tp 8 G S 1 . Then either 
s = 1 or s = 2 (see. Fig. 4 a to the Example 2). 



c = e 3x + e 2x+y +e x-y 



91=71/4 










\ / x 






\ 







cos(v??) = sin(v??) = ^= 
cos(y?!J) + 2sin(y^) = 
2cos(<p§) + sin(</?°) = 
J W}(y?) > 
Fig. 6 a. 



Zones of positivity for 
( x>0 

hwj}(v) ■ I y> -2x 

Fig. 6 b. 



We choose in the class {ce 1 ^} a representative d with minimal set of 
zeroes of the indicator of growth Wj — » + W = Wj for c' 



w 
e vv c. 
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Possibility 3 with 2 opposite zeroes Iw i'Pi) — Iw (v 9 + 7r ) i s realized for 
c = e Wl + e _H/l . It is impossible to reduce it to the case 1. If s = 1, then 
it is possible to reduce it to the "Stable Positive Case" case Iwj+wif) > 
by choosing W. The last case we call "stable". The function d in this 
class {Wj + W = Wj} where I w >(ipi) = has an isolated zero, we call "a 
boundary function". It is easy to see that in the stable case where the 
function c > grows exponentially in all directions, its opposite c~ x has 
exponential decay in all directions. 

Let d = Co G {ce 1 ^} be stable (i.e has exponential growth in all directions). 
All such functions c are the ground states for L_ where L p = L + ©L_. 
All stable functions c" = cw from the class {c e w } define the ground 
states for L p located in the sector L_. 

e -i(ay-Px)/2 
e -i(ay-f3x)/2 

y/cw 

Conclusion. The operator L_ has the family of the square integrable 
ground states ^w, parametrized by the linear forms W = ax + fly (or 
by the pair a, (3) G R 2 such that the set {Wj + W} has strictly positive 
indicator of growth Iwj+w > everywhere on the circle. This domain T 
on the plane R 2 with coordinates (a,/3) is a convex polygon T. The 
interior of T is nonempty for all stable cases (i.e. if the set of linear forms 
Wj — W s generates linear space of dimension more that one). Inside T 
the ground states belong to C 2 (M 2 ), W G Int T. On the boundary of 
T the states ^ w do not belong to £ 2 (K 2 )- Therefore they represent the 
bottom of continuous spectrum for the operators L_ , and L p = L + © L_ . 
Apparently, the lowest level for the another sector L + , is strictly positive. 

Every state for L_ except Co generate a nonzero current Jw defined 
by the phase of the complex function ipw The Current vectors Jw cover 
some convex bounded domain on the plane. 

Important fact. The magnetic flux is divergent on the plane R 2 : 

Bdxdy = ■ <j> I{ Wj y((p)d(p + o(R), R — > oo, 

x 2 +y 2 <R 2 




20 



where I{w} is the indicator of growth (see Appendix). 
II. For the generating functions of the form: 

c=J2 K 3 eWj + eWR ' q K cos W i,i + «9 sin ^,3) 
i 1 

we define the indicator of growth using only its part: the real subset 
should be chosen such that Kj > 0. It defines a domain T + 

as above. (Possibly, it is necessary to change the sign c — > — c, if this 
operation leads to the bigger domain). Anyway, all other linear forms 
Wj,WR jQ should belong to the domain T + defined by the main positive 
part of linear forms selected above. Under these conditions there exists a 
nonempty set of coefficients such that c ^ 0. It is possible to choose 
remaining coefficients such that zeroes of c appear. 

2.2 Solutions of genus g=l. Operators with Bohm- 
Aharonov Singularity and the Magnetic Bloch func- 
tions 

In the current section we study the case of elliptic curves (g = 1). As w are go- 
ing to show, in this case we are facing a new very interesting phenomena, con- 
necting out paper with the case of non-zero magnetic flux [31 H] (magnetic 
translations and topological phenomena are also discussed in [6j El [9| [TO]). 
Let 

r = r" = c/a, 

where 

A = {2mio>i + 2m 2 W2, mi, m 2 G Z} C C 

denotes a lattice. We assume that u\ = 1 and the lattice is invariant with 
respect to the complex conjugation 

A = A 

(it is true, in particular, if co> 2 e iR). 

Let 00 1 = G I"", oo 2 = G r". Assume, that Q , . . . ,Q n G T' and 
R , . . . , R n G r" correspond to the intersection points V fl T". 

Let us express the ^-function through the Weierstrass cr-function and 
^-function . 
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We know that ((w) is meromporphic in C with first order poles at the 
points of A and 

C(w + 2u s ) = CH + 2rj S: (1) 

where r] s = ((uj s ). 

The lattice is invariant with respect to the complex conjugation, therefore 



CW = CH- 

The function a(w) is analytic in C with first-order zeroes at the points of 
A and 

a(w + 2u s ) = -e 2ris{w+Ws) a{w), 

a(w - 2co s ) = -e- 2v ^ n) a(w). 

The invariance of the lattice with respect to the complex conjugation implies 
that 

a{w) = <j{w). 

The function ip" = tp\ T n has the following form (z,z£C,j)6 T", P = D")\ 

^ [P,) a(z + P)a(p + Py 

The function ip' = tp\ r i has the form 

nk, t> = e-«« ( ^ + ^^-Q,K p ^-Q„) ( 2) + . . . + 

V a(k + Pi) . . .a(k + P n+1 ) 
a(k + z + A n )a(k - gg) . . . o-(fc - Q w _i) \ 

(7(fc + P 1 )...(7(fc + P B+1 ) M*,*^, 

where 

^0 = Ql + • • • + Qn + Pi + • • • + Pn+1, 
A 1 = Q + Q 2 + --- + Qn + Pl + --- + Pn+1, 



A n = Q + --- + Q n _! + Pi + • • • + P n+1 , D' = P\ + h Pn+L 

The compatibility conditions 
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imply 



SC(Qo) a (Q° + Z + A)MQo ~ Qj) ■ ■ ■ a (Qo - Qn) f / 

a(Q + P 1 )...a(Q + P n+1 ) ' ' 



a(z + P)a(Ro + P)' 

-O A 

fn(z,z) 



-zC,{Q n ) a (Qn + z + A n )a(Q n - go) . . . a(Q n - Q n -i 



a(Q n + P 1 )...cr(Q n + P n+1 ) 

e -zC(R n ) °( R n + Z + P) 

a(z + P)a(Rn + P) ' 

therefore f s (z, z) : 

f = e -zC(R )+zC(Q ) a ( R + Z + P) „ 

a(z + Q + --- + Q n + P 1 + --- + P n+1 )a(z + P) °' 

( 7(g + p 1 )...cr(go + Pn+i) 

D 



a(P + i>(Qo - Ql) • • . ^(Qo - Qn) 



f = e -zC(Rn)+zC(Q n ) (P TO + Z + P) 

^ ( T(^ + g + --- + Qr i + Pl + --- + Pn + l)^(^ + P) 

a(Q n + P 1 )...a(Q n + P n+1 ) 



a(R n + P)a(Q n - Qo) . . . a(Q n - Q n _^ ' 
We have 

= /o( ,,,-) + ... + 

V 0"IPJ • • • ^IPi+lJ 

<t(z + A n )<r(-Q ) • • • o-(-Qn-i) , , _A 

Jn{Z,Z) . 



a(P 1 )...a(P k+1 ) 

Let us point out that all f s have the same factor at the denominator 
a' a" = a(z + Q + ■ ■ ■ + Q n + Pi + • • • + P n +i)cr(z + P). 

Multiplying all f s to a' a" we obtain nonsingular functions c = ca{z + Q + 
D')a(z + D"). The corresponding functions ip = ipa(z + Q + D')a(z + D") do 
not have Bloch-Floquet properties. They associated with the ground state of 
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Pauli operator with magnetic field B = 1/2 A (In c). We are going to discuss 
their analytic properties and relationship to the Dubrovin-Novikov states 
(1980) below. 

Easy to prove following statements. 

Proposition 1: All solutions c can be presented in the form 

c = ^2® q exp{-z((R q ) + zC(Q q )}o-(z + R q )a(z - Q q ) 
q 

Here a q , R q , Q q can be any generic set of numbers, no restrictions. 

Proposition 2: Every real solution can be presented as a sum of two type 
of terms like in Proposition 1: The Type 1 where a q G K and R q = —Q q , 
and The Type 2 where we have pair of indices j, I satisfying to relations 
ctj = on, Rj = -Qi, Ri = —Qj. 

Definition. We say that real solution c has Type (k, I) if it is 
realized as a sum of k Type 1 terms and I Type 2 terms (total 
number of intersection points is k + 21). 

No problem to choose parameters leading to the nonzero real function 
c^O.cGt. 

How to choose c with periodic magnetic field B with same pe- 
riods as our lattice? 

Example. Let us take n = 2 (3 intersection points) and real function c 
of the type (1, 1): 

c = aexp{zi] + zfj}\a(z — co)\ 2 + 

+{3exp{z((Q 1 ) - zC(Ri)}o-(z - Q 1 )a(z + Rj) + CC 

Here CC means "complex conjugate". We can choose parameters such that 
this expression is positive everywhere. 

Proposition 3. Fix a G R and (3 G C. There exists a countable number 
of data R2 = —R\ = —Qi = Q2,Ro = = —Qo such that magnetic 
field B = 1/2 A log c is periodic with the same periods as lattice. This set is 
determined by one number which satisfies to the equations 

U — U = —inn, V + V = —nm, m,n G Z 

17(A) = A77 — wC(A), V(X) = At/ - w'C(A) 

Sum of such positive expressions c = c s also leads to the Algebro- Geometric 
operator with nonsingular real magnetic field B with flux equal to one quan- 
tum unit. 
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Conjecture: Every periodic nonsingular magnetic field B with 
flux equal to one quantum unit through the elementary cell, can be 
approximated by the Algebro-Geometric Fields described above. 

Let us make comparison with the Dubrovin-Novikov bases [3, 4|. 

Our function ip is an eigenfunction for the Pauli Operator L p with the 
magnetic field B. Is it a magnetic analog of the Bloch functions? A complete 
family of the Magnetic Bloch Functions, parametrized by the points of the 
torus T 2 = I"" was constructed in [31 H] for all non-zero values of magnetic 
flux for an arbitrary nonsingular periodic field B(x,y). We consider at the 
moment the case of flux equal to one. The whole family determines a direct 
summand in the Hilbert space £2(K 2 ). These functions were written in the 
form: 

ip = (const) exp{S}cr(z + a') 

where — 1/2 AS = B, and inversion of the Laplacian was especially normal- 
ized. Here S = S + (const) z, where S is a real functions, independent on a' 
and (const) depends on a' and the lattice. In our case we start with function 
ip" holomorphic in z because the actual ground states (corresponding to the 
spectrum in the Hilbert Space) belong to the spin-sector described by the 
curve T": We have to reduce the operator L p = L + @ L~ to the self-adjoint 
form by the "Gauge Transformation". We choose the proper spin-sector 
where the ground states are located, i.e. eigenfunctions decay at infinity (as 
for genus zero). It corresponds to the curve T". For that sector we need to 
divide ip — ty" by the y/c. Use here notations a = P,b = A,d = p — P. So 
we have finally 

<j(z + d — a) 



1>(p,z) = (1/V~c)r(p,z) = (xJa(z-a)a(z-b))/(Vl)e-^ 



a(z — a)a(d 



= (const) a / a( f h \ (l/Vc\e-« d **<T(z + d-a) 
y a(z - a) 

where p G T" = C/A. In order to remove a singular Bohm-Aharonov 5-term 
we apply a unitary singular gauge transformation and multiply the result to 
a linear exponent 

V a ( z - 0) 
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where (const) depending on p is chosen to obtain bounded magnetic Bloch 
eigenfunctions. 

The last function if) leads exactly to the Dubrovin-Novikov magnetic 
Bloch family of eigenfunctions for the smooth field B: 

$ = (const) [l/yfc)e-^ d)z e {const)z a(z + d - a) = exp{S}a(z + a') 

For smooth fields this family is uniquely defined. 

Another form of this argument is following: we can always multiply ip" 
by any holomorphic function, which we choose as a (z - a )e( const > with cons 
depending on p. The product also satisfies to the Cauchy-Riemann (=" self- 
duality") equation as an "Instanton" in the first nonself-adjoint form of the 
Pauli operator. After that we multiply result by the 1/y/d realizing a non 
unitary "Gauge Transformation" leading to the self-adjoint form of L p in the 
sector where ground states are located. Finally we get exactly the Dubrovin- 
Novikov Magnetic Bloch function (see [3]). 

These arguments work because we are dealing with the " Instanton Fam- 
ily" (i.e. satisfying to the first order Cauchy-Riemann equation dip" = 0) for 
the operator L p written in the first non-selfadjoint form. It satisfies to the co- 
variant first order equation [8+ l/2d(lnc)][(l/\/c)ip"] = after nonunitary 
gauge transformation realized as a division of all eigenfunctions including 
ground-states by the factor y/5. Here a, 6, d are constants with the proper 
reality restrictions replacing the points a = P,b = A = J2 S — S s Qs,d = 
p — P correspondingly on the elliptic curve. In particular b — a. 

For g > 1 we always have magnetic flux with more than one 
quantum unit. We have to use Riemann surfaces V with selected 
point ooi such that solutions of the corresponding KP hierarchy 
is Elliptic in the variable x in order to get magnetic fields B and 
B periodic in both directions in the z-plane. The theory of elliptic 
solutions to the KdV equation was started by Dubrovin and Novikov in 
1974. For the KP hierarchy it was developed by Krichever since 1979 . A 
number of works were dedicated to it in the later literature. Details will be 
presented in the next work of the present authors. 

Let us discuss here an extremely important physical question: 
What is a quantum Bohm-Aharonov Phenomenon? How the 5- 
term in magnetic field affects the spectrum? 

In our case the singular magnetic field B with the 5-type singularity has 
an algebro-geometric realization. It has a zero magnetic flux. The family of 



26 



the complex Bloch-Floquet eigenfunctions is found for it which has very spe- 
cific analytical properties valid only in the case of zero total flux through the 
elementary cell. It was explicitly calculated in the appropriate spin-sector. 
Its "instanton part" is simply the Baker- Akhiezer function ip" . After reduc- 
tion to the self-adjoint form of operators L p it became ip" / ' yfc with same 
quasimomentum. Other part ip'/y/c certainly does not have the instanton 
form outside of the intersection points. Is everything correctly and uniquely 
defined for the singular operators of this type? Are the Bloch-Floquet mul- 
tiplicators (whose logarithms divided by periods define the components of 
quasimomentum) canonically well-defined for such singular operators? 

This question needs clarification. Our statement is following: This 
family is correctly defined as a limit of corresponding families for 
the smooth operators with zero magnetic flux. Such a procedure 
to define spectrum and the whole complex family of Bloch-Floquet functions 
can be realized for the zero level by the family of Riemann surfaces T T de- 
generating to r = r'P|r". As we can see, only some special isolated state 
in this family might serve the standard Hilbert Space. In the Hilbert space 
£2(K 2 ) such state corresponds to the bottom of the continuous spectrum for 
the field B. So either the spectrum near is continuous or the dispersion 
relation near the zero point is identically trivial, and we have in fact more 
Bloch functions on the level e = 0. It is exactly the case here. We will clarify 
this question below. 

The Bloch-Floquet multiplicators of the family ip"/ y/c £1X6 K>x •> ^"ii 

They 

are equal to 

k x = exp{— 2u((p) — 2i]p}, K y = exp{—2i((p)u' — 2ii]'p} 
The equations \k x \ = 1,\k v \ = 1 can be easily solved 

P R = -C(P)V?7,/ = -CV/eta' 

or 

— C(p) — p R i]/ UJ + ip Ir ri' /u' 

So they are non-unitary \k\ ^ 1 for other points p G V" . 

Removing singular part from B we are coming to the Magnetic- 
Bloch Functions found in 1980: the multiplicators (i.e. the eigen- 
values of the Magnetic Translations) became unimodular \k\ = 1. 
Remaining spectrum will be separated from by the finite gap for 
the field B. 
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The formal procedure is following: We introduce function 

^new(p,z) = (lp"/\/c)exp{u(p)z} 

choosing u(p) such that all multiplicators became unitary k — > k: 

I 2u(p)o)| _ -I _ | p2i«(p)o)'| 
\' v x^ I — — \ y ! 

For the new multiplicators (after this renormalization of eigenf unction) we 
have 

k x = exp{2ip I [r) > u> — 770;'] /u/}, & y = exp{2ip R [?7u/ — 

Finally we remove singularity by multiplication ip new — > Ic'lVw™ replacing 
c by the smooth c in the formula above. So we are coming exactly to the 
magnetic Bloch eigenfunctions of Dubrovin and Novikov(1980) with multi- 
plicators k x , K y (may be shifted by constant which is inessential). 

The multiplicators k form together a point of the 2-torus k G T 2 forming 
the whole component space of the real quasimomentum at the zero level and 
nearby. We have a big complex 2-dimensional manifold M 2 of the Bloch- 
Floquet eigenfunctions for the singular operator L p with zero flux of the 
form 

m 2 = r" x cp 1 

compactified at the infinities u G CP 1 = C|Joo. It presents an irreducible 
component of the whole Bloch-Floquet manifold for this operator. The dis- 
persion relation e — > C degenerates at this component, i.e. e = identically. 
The Bloch eigenfunction already written above has a form 

ty(p, u, z) = ip" ' /y/c x exp{uz} 

where (p,u) G M 2 . This manifold presents exactly one component of the 
limit of the whole Bloch-Floquet manifolds M 2 for r ->■ and M 2 = M 2 [j M'. 
Every small purely magnetic perturbation leading to the smooth magnetic 
field B T closed to our singular field B, has quite similar Bloch-Floquet eigen- 
function: The Bloch manifold here is xCP 1 with the instanton-type Bloch 
eigenfunction like \I/ T = ip"(p, z)j \fc^ x exp{uz} where the curve might 
have an infinite genus. So the small electric perturbation of operator L p is 
needed. We perturb by the small periodic potential tU(x,u). Already the 
first order in coupling parameter r probably leads to a nontrivial dispersion 
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relation as a complex meromorphis function on the same manifold. It has 
the order r, so we have 

e T : M 2 ->■ C 

The real levels e T = const e it! give a function on the real " quasimomentum" 
torus T 2 whose levels are the real Fermi-curves. Its minimum lies nearby of 
the initial point on the curve T" found above where k x , k v are imaginary. 

We are going to calculate this perturbation in the next work. The complex 
level curves should have analytical properties at infinity typical for the Baker- 
Akhiezer functions ( maybe of the infinite genus where small handles appear 
from the "resonance points"). 

So it looks like the delta-term does not affect deeply the spectrum near 
the ground state. 



3 Appendix: The Asymptotic of Magnetic 
Flux 

Let us calculate the Asymptotic of Magnetic Flux through the round ball of 
radius R for the purely exponential case 

e 2 * = c = Kje Wj ,Kj > 0, Wj = R{aj cos(0) + sin(0)) 

j 

We have for the magnetic field B = — (A(lnc))/2. For the vector- 
potential restricted on the circle r = R in polar coordinates, we obtain 

A = § y dx - § x dy = ~R\5~] Kje Wj (aj cos(0) + bj sm(<f>))]d<f>/c. 

3 

Our assumption is that there exist exactly iV 

indices j = 1,2, ...,N such that the Indicator of our family I^ w ^(<p) = 
maxjlw {4>) where Iy/ = max[aj cos(0) + /3j sin(0), 0] is strictly positive, and 
all other indices p ^ 1,2, ...N are inessential (i.e. corresponding linear forms 
W p are located strictly inside of the convex domain T C R 2 with coordinates 
a, (3 numerating the rapidly decreasing ground state vectors of our operator). 
There are domains Aj on the circle S 1 , where I{w k } = Iwj{<j>) with end points 
Aj = [0°, 0jLi], and for j — N, 1 we have $ = (/)%. So A N is a neighbor of 
Ajv_i and Ai (i.e. our numeration ic circle contr-clockwise). 
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Our claim is following: Following Asymptotic Formula is true: 



N 



= J2 R ~ S J2 X ?^ S ^ + ^K 1 )(" 1 ) S } + (Remainder). 

S>1 j = l 

Apparently, this series is nonconvergent, since the coefficients Q s grow, as 
we think, as s\. About the remainder we claim now that its decay is more 
rapid than any negative degree of R. We claim only that the "Regularized 
Flux" 

I I Bdxdy+ l -R<l I {Wj} {4>)d4> = O (i 

is tend to zero in this sum, for R — > oo. Performing this calculation 
near the critical points 0°, we use following functions 

(W j+1 - Wj)/R = (a j+1 - a k ) cos(0) + (0 j+1 - pj) sin(0) = tj(z). 

Here = (0° + z), \z\ < e. It is located near the points 0° or z = 0: in this 
point Wj = Wj+i,tj = 0,z = 0, and the inverse function z(tj)) is given by 
the inverse series with a finite radius: 



z = 

k>l 



k>l 

d<f) = dz = J2 X f )tk 3 dt 3 



k>0 

We define numbers 



POO 

Q k (a) = / [aw k e' w /(l + ae- w )]dw 
Jo 

useful for the investigation of the difference 

/ a + \rI I {Wq} {<t>)d<t>. 
Js^ 1 Js 1 

Probably, Q k ~ k\. Our function c has exponential growth everywhere, but 
magnetic field has decay only outside of the small domains surrounding the 
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"critical" points 0^. It is easy to see that our vector-potential A after ex- 
tracting the Indicator of Growth i2i{vy }(0)(<i0), became exponentially small 
outside these small domains. Only two exponential terms 
in c are essential in every such small domain, between Aj and Aj + ±. Drop- 
ping all other terms in the sum for c = J2 q K q e Wq and for A in every such 
small area costs us exponentially small. In the area G A q we multiply 
both-numerator and denominator in the expression for ^4-by the exponent 
K~ l e~ Wq . The exponent e~ W] is the vertex of the convex polygon T contain- 
ing all functions d = ce w G T such that {ce Wq )~ 1 / 2 are the ground states of 
the Pauli Operator. We need q = j for < 0° (or G Aj), and q — j + 1 
for G Aj + i. So only two terms remain in the numerator and denominator. 

Similar result we obtain in the domain Aj just below the point 0° with 
inverse constant Kj + \/ Kj and exponent e^ Wj+1 ~ Wj ^ , plus we have to turn back 
the direction of integration. Taking e such that Re = 0(R S ),5 > 0, we see 
following: The integration between the local limits [0° — e, 0° + e] of such 
expressions with w = Rtj, which appear in our calculation of the regularized 
magnetic flux, can be extended to the limits [— oo,+oo]. It is true because 
the remaining terms have order 0(e~ R ): more precisely their decay is more 
rapid than any polynomial. 

Expressing the variable z = — 0° by the variable tj = (Wj — Wj +1 )/R, 
we are easily coming to our result. In the final integration we have a sum of 
integrals looking like 

POD 

Q s (a) = R' 3 - 1 / ae~ w /(l + ae- w )w s dw, 
Jo 

where w = ±Rtj. The sign is + and a = Kj+i/Kj for z < 0, and sign — and 
a replaced by gT 1 for z > 0. So we are coming to our result. 
Note, it is easy to show that the expressions 

POD 

Q' k = / e~ w w k dw 
Jo 

grow as (k\). Probably, it is true for our expressions Qk- 
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